Spectral functions, Fermi surface and pseudo gap in the t-J model 



P. Prelovsek and A. Ramsak 
Faculty of Mathematics and Physics, University of Ljubljana, 1111 Ljubljana, Slovenia 
J. Stefan Institute, University of Ljubljana, 1111 Ljubljana, Slovenia 
(September 4, 2001) 

Spectral functions within the generalized t-J model as relevant to cuprates are analyzed using the 
method of equations of motion for projected fermion operators. In the evaluation of the self energy 
the decoupling of spin and single-particle fluctuations is performed. It is shown that in an undoped 
antiferromagnet (AFM) the method reproduces the selfconsistent Born approximation. For finite 
doping with short range AFM order the approximation evolves into a paramagnon contribution which 
retains large incoherent contribution in the hole part of the spectral function as well as the hole- 
pocket-like Fermi surface at low doping. On the other hand, the contribution of (longitudinal) spin 
fluctuations, with the coupling mostly determined predominantly by J and next-neighbor hopping 
t' , is essential for the emergence of the pseudogap. The latter shows at low doping in the effective 
truncation of the large Fermi surface, reduced electron density of states and at the same time 
quasiparticle density of states at the Fermi level. 



PACS numbers: 71.27.+a, 72.15.-v, 71.10.Fd 
I. INTRODUCTION 

One of the central issues in the experimental and 
theoretical investigations of superconducting cuprates is 
the understanding of low-energy electronic excitations 
in these compounds Q, where the clue to their high- 
temperature superconductivity lies possibly already in 
their anomalous normal-state properties. In recent years, 
in particular, the remarkable progress in the angle- 
resolved photoemission spectroscopy (ARPES) experi- 
ments revealed quite universal development of electron 
spectral properties as a function of doping. In most in- 
vestigated Bi 2 Sr 2 CaCu 2 02 +< 5 (BSCCO) ARPES shows 
quite well defined large Fermi surface in the overdoped 
and optimally doped samples at T > T Cl whereby the 
low-energy behavior with increasing doping in the over- 
doped regime qualitatively approaches (but does not in 
fact reach) that of the normal Fermi-liquid with under- 
damped quasiparticle (QP) excitations. On the other 
hand, in the underdoped materials the QP dispersing 
through the Fermi surface (FS) are resolved by ARPES 
in BSCCO only in parts of the large FS, in particular 
along the nodal (0, 0)-(7r, tt) direction J2|, indicating that 
the rest of the large FS is truncated j|, i.e. either fully or 
effectively gaped. At the same time near the (tt, 0) mo- 
mentum ARPES reveals a hump at ~ 100 meV [Q, which 
indicates on the existence of the pseudogap scale, which 
is consistent with the characteristic temperature pseu- 
dogap scale T* > T c , which appears also as a crossover 
in several other quantities: uniform susceptibility x(^)> 
resistivity p(T), the specific heat Cy(T) and the Hall 
constant Rh{T) Q. Although the latter anomalies in 
thermodynamic and transport quantities are quite simi- 
lar (or even better confirmed and more pronounced) also 
in La2- a: Sr a ;Cu204 (LSCO), spectral properties of the 
latter are qualitatively different from BSCCO, pre- 



sumably due to the crucial role of stripe structures in 
the LSCO in the regime of intermediate doping. 

Since electron spectral functions are the central quan- 
tity to characterize and understand anomalous Fermi liq- 
uid in cuprates, they have been the subject of numerous 
theoretical studies. There appears to be at least a theo- 
retical consensus on the spectral functions in an undoped 
reference antiferromagnet (AFM), describing a single 
hole inserted in an AFM behaving in two-dimensional 
(2D) planar system as a QP with strongly renormal- 
ized mass and large incoherent component. The spectral 
function is well captured within the selfconsistent Born 
approximation (SCBA) || for the simplest relevant t-J 
model, whereby for an agreement with experiments on 
undoped cuprates, i.e., Sr 2 Cu02Cl2 || and Ca2Cu02Cl2 
longer range hopping terms have been invoked 

For larger (finite) doping in the lack of reliable ana- 
lytical techniques numerical approaches have been used 
extensively. Starting with the simplest models for corre- 
lated electrons in cuprates, i.e., the Hubbard model and 
the t-J models on planar lattices, numerical studies em- 
ploying mainly the exact-diagonalization and the quan- 
tum Monte Carlo methods were able to confirm some 
gross features consistent with experiments. In particular 
this includes: a) the existence of the large FS already in 
moderately doped AFM |p^|,p[, b) the overdamped char- 
acter of QP at the intermediate doping J]~T|,|T^ ] , consistent 
with the marginal Fermi-liquid concept pl3|, c) the pseu- 
dogap features at lower doping in spectral functions |l4| ] 
and in the density of states (DOS) jl2|], and d) quite vis- 
ible contribution of longer range hopping Whereas 
numerical studies confirm the relevance of simplest mod- 
els and results for cuprates, still they are prevented to 
approach effectively the low-energy regime and results 
need a proper analytical (phenomenological) interpreta- 
tion. 



1 



Analytical approximations to spectral properties of rel- 
evant models for 2D systems at finite doping have proven 
to be very delicate due to strong interaction between elec- 
tron excitations, spin degrees and pairing fluctuations. 
For the one-band Hubbard model spectral functions have 
been evaluated within the random-phase approximation 
for AFM fluctuations |l5| and within the self-consistent 
conserving (FLEX) theory jl6| , both restricted to modest 
U/t. Starting with the t-J model strong correlations are 
explicitly taken into account in slave boson theories |l7| , 
where it is however difficult to incorporate AFM spin fluc- 
tuations. The latter play the essential role in phenomeno- 
logical theory of the spin-fermion model leading to the 
nearly AFM Fermi liquid fig}] . Recently aspects of the 
pseudogap features in the underdoped regime have been 
found in this model evaluating the self energy involving 
a strong coupling to AFM spin fluctuations fl9|j2C| ] . 

Concerning the origin and the explanation of the pseu- 
dogap scale is seems plausible that at low doping it is 
related to the exchange J since T* ~ J in low-doping 
materials, whereas it merges T* ~ T c in optimally doped 
samples. This indicates on the importance of AFM spin 
correlations for the emergence of the (large) pseudogap 
as found also in the numerical studies JT4| , pT|JT^ ] and in 
phenomenological model studies [ p^|j20[ . The rcnormal- 
ization group studies of the Hubbard model |pl[ (with 
moderate U/t) also reveal the instability of the normal 
Fermi liquid close to the half-filled band (insulator) and 
a possible truncation of the Fermi surface, but here ei- 
ther the spin or pairing fluctuations can take the dom- 
inant role depending on FS nesting conditions whereby 
the longer-range hopping terms can be quite important. 

One of present authors introduced the equations-of- 
motion (EQM) method (2^] for the evaluation of the spec- 
tral functions within the t-J model. It has been shown 
that EQM for projected fermionic operators implicitly 
reveal an effective spin-fermion coupling. Using the sim- 
plest decoupling it was possible to relate the overdanxped 
marginal-type character of QP to the marginal dynam- 
ics of spins [^2| but also to treat the superconducting 
fluctuations and transition j23). EQM method has been 
also applied for the Hubbard model [Q. The analysis 
of spectral functions within the t-J model has been re- 
cently improved [ p5| by more appropriate treatment of 
the self energy by dealing separately with: a) the strong 
coupling to transverse short-range AFM spin fluctuations 
- paramagnons and b) the moderate coupling to slower 
longitudinal fluctuations of the AFM order parameter. 
The theory and results were partially presented already 
in Ref . |2^] . The aim of this paper is to present the theory 
and in particular its consequences and results in more de- 
tail. The main advantage of this theory is that it allows 
for the study of the evolution of spectral function as a 
function of doping: from an undoped system up to mod- 
erate doping. The emphasis is on the results at T = 0, 
in particular on: a) the development of the FS from a 



hole-pocket like into a large one, b) the emergence of the 
pseudogap in the spectral function and related effective 
truncation of the FS most pronounced near (n, 0) in the 
Brillouin zone, c) anomalous properties of QP at the FS 
remaining well defined even in the pseudogap regime, d) 
depleted DOS and moreover the quasiparticle DOS (re- 
lated to specific heat coefficient) with doping. 

The paper is organized as follows: In Sec. II EQM 
method for the spectral function within the t-J model 
is summarized. Sec. Ill is devoted to the evaluation of 
the self energy within the decoupling approximation sep- 
arately treating the paramagnon contribution E pm and 
the contribution of longitudinal spin fluctuations Ejf . In 
Sec. IV results of the simplified analysis of the pseudo- 
gap features are presented taking into account an effec- 
tive renormalized band and explicitly Eif . Sec. V presents 
results of the full selfconsistent solution for spectral func- 
tion as a function of doping. 

II. EQUATIONS OF MOTION 

In order to take explicitly into account strong correla- 
tions we study the t-J model, 

H = -^2 hjc] s Cis + J ^(Si ■ Sj - ^rurij), (1) 

i,j>s (ij) 

where fermionic operators are projected ones not allowing 
for the double occupancy of sites, i.e., 

4 = (i-^,- s )4- ( 2 ) 

Since longer range hopping appears to be important for 
the proper description of spectral function in cuprates, 
both for the shape of the FS at optimum doping materials 
H as well as for the explanation of ARPES of undoped 
insulators we consider besides tij = t for the n.n. 

hopping also ty = t' for the n.n.n. hopping on a square 
lattice. 

Our goal is to evaluate the electron Green's function 
(propagator) directly for projected fermionic operators, 

G(k,w) = ((ck a ;4 a ))w = 

= -i e^+^({c ks (t),cl} + )dt, (3) 
Jo 

which is equivalent to the usual propagator within the 
allowed basis states of the model, Eq. ([!]). In the EQM 
method p^j one uses relations for general correlation 
functions 

u((A;B)) u = ({A,B} + ) + (([A,H];B)) U = 

= {{A,B} + )-{{A;[B,H])) UJ . (4) 

and applies the propagator G(oj) = ((A;A'))u- If we de- 
fine the (orthogonal) operator C as 
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[A,H] = QA-iC, {{C,Ai}+) = 0, (5) 
we can express 

G(u) = G q (uj) + ^G ( W ) 2 ((C;Ct)) w , 



Go(w) 



(6) 



Identifying the self energy E(w) as the irreducible part 
of ((C; C^)) w we can express Eq. (||) as 

GH = 7~~yT7\' S(o;)~i<(C;Ct))*r-. (7) 

Within the diagrammatic technique E(w) corresponds to 
the contribution of irreducible diagrams. Generally E(w) 
can be defined as a memory function within the Mori pro- 
jection method ^7|. Anyway in most cases the successful 
application of the method relies on the appropriate de- 
coupling or other approximation of the memory function 
E(w) @. 

Applying the formalism to the propagator, Eq. (j^), 
we have to deal with the EQM for c ks with a nontrivial 
normalization factor, 

a = jj E({ £ -> 4i+> = i - f = ^(i + ch). (8) 

By taking explicitly into account the projection in Eq. (|^) 
the EQM follow, 

[Cis,H}= - ^ hi K 1 ~ n i,s)cjs + S^Cj-s] 

i 

+ Jj ( 2sS ]C™ + 25jci,_ s - nj-Ci.) (9) 

j n.n.i 

with s = ±1. We express 'bosonic' variables in to terms 
of spin and density operators, i.e., n^_ s = rij/2 — sSf . 
Assuming as well that we are dealing with a paramag- 
netic metal with (Si) = and a homogeneous electron 
density (rii) = c e , we get 



*Cks — [C ks , -ff] — CkCks 



(10) 



and 



[c ks ,ff] = [(l-|) £k -J Ce ]c ks + -^£ 



,"i kq 



1- - 



[s5qCk_q, s + S^Ck-q,-* - -n q C k _q jS ] , (11) 

where = rif — c e , m kq is the effective spin-fermion 
coupling, 



m-kq = 2J7 q + e k _ 



(12) 



and e k is the bare band dispersion, i.e., for the model (Q) 
on a square lattice 



^ = -4t 7 k - 4i' 7k , 
7 k = — (cos k x + cos k y ), 7 k = cos k x cos fc y . 



(13) 



Equations (||) , (|Tl|) , (|l3|) define also the 'renormalized' 
band 



Ck 



-<{[c ks ,ff] ) c ks }+) 
C - 4?7it7k - 4?72i'7 k , 
a H (S ■ Sj), 



(14) 



where 77^ are determined solely by short range spin cor- 
relations and £ is a k independent term (still dependent 
on various static correlations). 

Above quantities determine the propagator 



u + A* - Ck - S(k,w) ' 



(15) 



and the corresponding spectral function A(k, to) = 
— (l/7r)ImG(k, u>), provided that we find a method to 
evaluate E(k, lu). 



III. SELF ENERGY 

A. Undoped antiferromagnet 

It is desirable that in the case of an undoped AFM 
our treatment of E and the spectral function reproduces 
quite successful SCBA equations || for the Green's func- 
tion of a hole in an AFM. Let us concentrate here on the 
relevant n.n. hopping, since the t' term represents a hop- 
ping on the same sublattice within an ordered AFM and 
is therefore nearly free. For the SCBA the reference state 
is the Neel state with rii S = 0, 1 for i = A, B sublattices, 
respectively, and the SCBA effective Hamiltonian can be 
written as 



H h = -tJ2( h * h ] a j + h jh\a\) + Hj, 

(ij) 



(16) 



where hi represent holon operators and cii spin flip oper- 
ators. The corresponding holon EQM then follows from 
Eq. © 

- ijh\ = [h u H h ] = t J2 h]{a\ + a 3 ). (17) 

j n.n.i 

It is now straightforward to establish the relation of 
Eq. ([IT]) with the EQM for Ci S by considering one Neel 
sublattice i = A with the reference state n.i S = 1. In this 
case 1— rii.- s = 1 and by formally replacing Cj s = c.j.- s Sj 
we get by considering only the t term in Eq. (Bh, 



ii—rC; 
dt 



t J2 (S? + Sf)c jt _ s 



(18) 
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To be consistent with the SCBA we neglect here the J 
term in Eq. (|9|) since anyhow J -C t. Within the lin- 
earized magnon theory EQM jT^),([T8|) are formally iden- 
tical, so we can furtheron follow the procedure of the 
evaluation of EAFM(k, lo) within the SCBA to reproduce 
in the first place spectral properties of an undoped AFM. 
In this case we do not try to improve the SCBA, since 
the latter approximation is simple and yields both qual- 
itatively and quantitatively good results consistent with 
numerical studies and experiments. For an ordered 2D 
AFM where relevant spin excitations are magnons with 
dispersion w q we therefore get 

E AFM (k, lo) = 1 J2 M^G(k -q,co + Wq ), (19) 



Mkq = 4t(w q 7k-q + Vq7k), 



with 



<2J 



2lo c 



-7q1 



'2 J 



2U! 



^q = 2J A /l- 7 2. 



(20) 



Since in a Neel state we have rji = and hence the 
renormalized band vanishes, i.e., £k = 0, we reproduce 
the usual SCBA equations for the hole spectral function 
within the t-J model. The inclusion of the n.n.n. hop- 
ping t' is also simple within the SCBA since within the 
Neel state it does not induce a coupling to spin flips in 
Eq. (|^) and therefore enters G(k, lo), Eq. (|l5|), only via 
the band term (k ~ £ — 4t'7^. It should also be noted 
that in contrast to the usual SCBA our procedure deals 
directly with the electron propagator and not with an un- 
physical holon one. Moreover it allows a straightforward 
generalization to the case of finite doping. 



B. Short range transverse spin fluctuations 

For finite doping ch > we assume that spin fluctu- 
ations remain dominant at the AFM wavevector Q = 
(-7T, 7r) with the characteristic inverse AFM correlation 
length k = 1/£,afm- The latter seems to be the case 
for BSCCO as well as YB 2 Cu 3 6+2; , but not for LSCO 
with pronounced stripe and spin-density structures with 
Q.SDW 7^ Q- For the former case one can divide the spin 
fluctuations into two regimes with respect to q = q— Q: 

a) For q > k spin fluctuations are paramagnons, i.e., 
they are propagating like magnons and are transverse 
to the local AFM short-range spin ordering. Hence it 
make sense to use Eqs.(|l9|),(|2"o|) to represent paramagnon 
contribution to the self energy restricting the sum to the 
regime q > n. 

b) For q < k spin fluctuations are essentially not prop- 
agating modes but critically overdamped so deviations 



from the long range order are essential. Alternative ap- 
proximations to E(k, lo) have to be used here, as dis- 
cussed in Sec. III.C. 

We should also take into account that the SCBA for- 
malism has been derived for an undoped AFM, i.e., for a 
hole spectral function at lo < where only (added) holes 
participate. Since we are dealing with ct > we take 
into account the scattering of hole-like {lo < 0) QP by 
replacing full propagator G in Eq. ( |l~9| ) by the hole part 



G = F(k,«) = ± f ■ 



du'A(k, 



LO — LO 



(21) 



However it is easy to see that an analogous contribu- 
tion should arise from the electron-like QP with lo > 0. 
At finite doping case we therefore generalize (at T = 0) 
Eq. ( |l9] ) into the paramagnon contribution, 



E pm (k, w) = i^ [M^G~ (k - q, lo + w q ) + 



q,q>K 



M"k+ q , q G + (k + q,o;-a; q )], 



(22) 



which would emerge from Eq. (|i8|). The consequence of 
Eq. (||^) is that in general Im£ pm (k, uj > 0) ^ so that 
also electron-like QP can be damped due to magnon pro- 
cesses. We do not consider here effects of T > which 
could be easily incorporated through the magnon occu- 
pation, but in most cases do not have a strong influence 
at low T < J. 

Here we stress two features of our approximation for 
paramagnon contribution S pm : 

a) we are dealing with a strong coupling theory due to 
t > ujq and a selfconsistent calculation of S pm is required, 

b) resulting S pm (k, lo) as well A(k, lo) as are at low dop- 
ing quite asymmetric with respect to lo — 0. As in 
an undoped AFM the hole part G~ with the weight 
oc (1 — Cfj)/2 ~ 1/2 generates a large incoherent part 
in A(k, lo -c 0). On the other hand, G + has less weight 
cx Ch and consequently the scattering of electron QP is 
in general much less effective. 



C. Coupling to longitudinal spin fluctuations 

Discussing the self energy at finite doping, Eq. j22] ) 
represents only one contribution and we have to recon- 
sider EQM We note that at c h > C ks con- 
tains a remainder of a 'free' term cx Ck s , which should be 
however neglected when evaluating the 'irreducible' part 
entering E, Eq. (0). Considering within the simplest ap- 
proximation only the mode-coupling terms in Eq. (^) we 
also neglect the coupling to density fluctuation n q which 
should contribute much less to E in the absence of charge 
ordering or charge instabilities at low doping. 
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Taking into account only spin fluctuations, at Ch > 
we are dealing with a paramagnet without an AFM 
long-range order and besides the paramagnon excitations 
also the coupling to longitudinal spin fluctuations become 
crucial. The latter restore the spin rotation symmetry in 
a paramagnet and EQM ( |TT| ) naturally introduce such 
a spin-symmetric coupling. Assuming within a simplest 
approximation that the dynamics of fcrmions and spins 
is independent, 

(55(t)p k - q) .(t)S« q 4_ q/)a ) ~ 
<5qq'(5 q (t)^ q )(c k _ q>s (t)cL qiS ), (23) 

we get for the contribution of longitudinal fluctuations 



£ lf (k,a,) = 5£ 



du>idu>2 



g{ui, w 2 ) 



A(k-q,o;i)x // (q, t ^) 



(24) 



2 L 2 



cth 



/3u> 2 - 



where \ is the dynamical spin susceptibility 



x(q»w) = -i 



[S'{t),S'])dt. 



(25) 



Such an approximation for £ has been introduced within 
the t-J model in Ref. p^ ]. However, quite analogous 
treatment has been employed previously in the Hubbard 
model flq] a nd more recently within the spin-fermion 
model yjf(]]. 

Several comments are in order to define quantities en- 
tering Eq. (p5| ): 

a) EQM (O) induce an effective spin-fermion coupling, 
which would emerge also from a phcnomcnological spin- 
fermion Hamiltonian with the coupling parameter mkq, 
Eq. (|l2|). In order that such a Hamiltonian is hermitian, 
the coupling should satisfy 



"•k-q,-c 



(26) 



which is in general not the case with the form Eq. ( |12[ ) 
therefore we use furtheron instead the symmetrized cou- 
pling 



m kq = 2J 7q + -(e k _ 



(27) 



Here we should point out that in contrast to previous re- 
lated studies of phcnomcnological spin-fermion coupling 
p5| , pj|j20| , our mkq (as well as mt q ) is strongly depen- 
dent on both q and k. It is essential that in the most 
sensitive parts of the FS, i.e., along the AFM zone bound- 
ary ('hot' spots) where k = |Q — k|, the coupling is in 
fact quite modest determined solely by J and t' . 
b) Since we are dealing with the paramagnetic state, all 
quantities should be spin invariant, i.e., the x Q ^(q, w) = 



Sapxfav). Since EQM (jLlh are invariant under spin 
rotations we have besides the S z term analogous terms 
with S~ , S + . Still we expect r s = 1 instead of r s = 3 
since only the coupling to longitudinal (to local Neel spin 
order) spin fluctuations is considered here, while the cou- 
pling to transverse fluctuations has been already taken 
into account by E pm . 

c) In Eif only the part corresponding to irreducible dia- 
grams should enter, so there are restrictions on the proper 
decoupling. We will be interested mostly in the situation 
with a pronounced AFM short-range order where longi- 
tudinal fluctuations are slow, i.e., with the characteris- 
tic frequencies lu k < 2Jk -C J. The regime is close to 
that of quasistatic xiQi 10 ) where the simplest and also 
quite satisfactory approximation is to insert for A the 
unrenormalized A , the latter corresponding in our case 
to the spectral function without Eif but with E = E pm . 
Such an approximation has been introduced in the the- 
ory of a pseudogap in CDW systems p9[] , used also in 
related works analyzing the role spin fluctuations [15|, 
]l9|] , and recently more extensively examined in Ref. [ 30 1 . 
In the opposite case of full self-consistent treatment with 
A = A we would overcount the influence of fluctuations, 
although the results would probably appear not so much 
different as shown on simpler systems ]30||. 



For %(q, w), Eq. (|25[), at q, > and possibly T > 
we do not have a corresponding theory, so we treat it as 
an input. On one hand, x(Q:^) is restricted by the sum 
rule, 



rcth(^) X "(q,u;)^= \(l-c h ) 
„ J ~ 



(28) 



At the same time, the system is close to the AFM insta- 
bility, so we assume spin fluctuations of the overdampcd 
form H 



X"(q,w) oc 



(q 2 + K 2 )(uj 2 +uj 2 K )' 



(29) 



The appearance of the pseudogap and the form of the 
FS are not strongly sensitive to the particular form of 
x"(q, oj) (at given characteristic k and w K ) provided that 
x(q, is not singular as, e.g., is the case of marginal 
Fermi liquid scenario [Q. It has been shown |^] that 
the latter form is needed to get generally overdamped 
QP with vanishing QP weight (at T = 0) in spectral 
functions at intermediate doping. 



IV. PSEUDOGAP ANALYSIS 

Full calculation of the spectral functions A(k, uj) within 
the presented theory requires a selfconsistent solution for 
E = E pm + Eif, where besides the model parameters 
t, t' , J and the doping c^, an input are also fi, k, rji, 772- 
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k, rji , r]2 are given by short-range spin correlations de- 
pendent mostly on c/j and can be taken from various an- 
alytical |3l| and numerical |32|,|| calculations within the 
t-J model. At the same time in a selfconsistent theory 
should be fixed via the DOS 



as 



c h = 1 - / f(u>)Af(iv)du 



(30) 



(31) 



Results of such selfconsistent calculation are presented in 
Sec. V. 

On the other hand, to establish characteristic features 
of the pseudogap and the development of the FS we follow 
first a simplified analysis. We namely notice that the 
effect of S pm are threefold: 

a) to induce a large incoherent component in the spectral 
functions at u> <C in particular at low and intermediate 
doping, 

b) to renormalize the effective QP band relevant to the 
behavior at u> ~ and at the FS, and 

c) to cause a transition of a large FS into a small hole- 
pocket- like FS at Ch < c* h <C 1, 

Result b) can serve as a starting point for the discus- 
sion of the pseudogap and FS features at finite doping. 
If we define the effective band as 



Zef 
k 



k -^[Ck + E pm (k,0) 
9E pm (k, u>) 



1 



j=0- 



we get for the effective spectral function 



(32) 



(33) 



which can be used to evaluate Sif. We restrict ourselves 
here to the regime of intermediate (not too small) doping, 
where e£ f defines the large FS. 

Let us concentrate on results for T = 0. The simplest 
situation where £if can be evaluated analytically is the 
quasi-static and single-mode approximation (QSA) which 
is meaningful if oj k <^ t, n <C 1. In this case we insert into 
Eq. (§) 

~ X "(q, w) ~ ~<5(q - Q)[8(w - v) - S(lo + v)\, (34) 
with v 
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0. We get 



Ejf^k.u,) 



4 u-4-<* 



(35) 



and 



Cfl SA %u) = 



aZ«(u> 



^ ) 



("-4f-q)("-3f)-A2 



(36) 



l kQi 



The spectral functions show in this approximation two 
branches of E^ , separated by the gap which opens along 
the AFM zone boundary k = Ilafm where e^-Q = e k • 
Since 7k^FAf = the relevant (pseudo)gap scale is 



^-^T S \2J -U'cos 2 k x \ 



(37) 



^k does not depend on 



It is instructive to realize that 
t, but rather on smaller J and in particular t' . For t' < 
the gap is largest at (jr, 0), consistent with experiments. 
Whether the (pseudo)gap appears at the Fermi energy 
uj = depends, however, on properties of e^ f . We 
do not expect that the gap opens along the whole AFM 
zone boundary, since in most cases e^ f crosses zero 
along (7r/2, 7r/2)-(7r, 0) so that within the QSA forms 
a hole-pocket-like FS. In fact, the results of the QSA are 
equivalent to the system with long range spin-density- 
wave order (AFM), where the doubling of the unit cell 
appears. 

Within the simplified effective band approach, 
Eq. (^), it is not difficult to evaluate numerically Ejf 
beyond the QSA, by taking explicitly x(q, uj), Eq. (p£|), 
for k > and io K — 2Jk. Integrals in Eq. ( p5| ) can be per- 
formed mostly analytically if we linearize the dependence 
of e£ within the relevant interval Sk < k. 

Let us for illustration present in this Section results 
characteristic for the development of spectral functions 
with most sensitive parameters k and [i, which both sim- 
ulate the variation with doping. We fix furtheron the 
model parameter J ft — 0.3 as relevant for cuprates. We 
take here t'/t = —0.3 close to values quoted for BSCCO. 
For simplicity we assume first that the effective band ejf 
is just renormalized e k (justified for an intermediate dop- 
ing, see Sec. V) with fixed values t e f/t = 0.3, t' ef /t = —0.1 
and Z cf — 0.4. More realistic treatment would require 
the variation of latter parameters with Ch but results re- 
main qualitatively similar. For k we take in accord with 
experiments JjJ and numerical results on the t-J model 

The choice of /i is somewhat more arbitrary since 
within an effective band approach the sum rule, Eq. (|3^) , 
cannot be used as a criterion. Nevertheless it is evident 
that /x determines the shape and the volume of the FS. 
In the following examples we choose fi such that at given 
k the DOS at the Fermi energy, Af(0), reaches a local 
minimum. This means that effectively the states near 
(71" , 0) are in the pseudogap and that the truncation of 
the FS is most pronounced (at given k). Such a choice 
of fi in fact also yields the volume of the FS (except at 
extreme c«l) quite close to the one consistent with the 
Luttinger theorem |33|. 
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In Fig. 1 we first present results for A(k, u = 0) at 
T = for a broad range of k = 0.01 — 0.6. Curves 
(evaluated at small additional smearing e = 0.02i) in 
fact display the effective FS determined by the condition 
G~ 1 (k.F,0) — 0. At the same time, intensities A(k,uj — 
0) correspond to the renormalization factor Zp. We can 
comment the development as follows. At extremely small 
k = 0.01 we see the hole-pocket FS which follows from 
the QSA in Eq. (36). In spite of small k the 'shadow' 
side of the hole pocket has smaller Zp. Already small 
k ~ 0.05 destroys the 'shadow' side of the pocket, i.e., 
the solution G _1 = on the latter side disappears since 
the singularity in Sif, Eq. (|35|), is smeared out by finite k. 
On the other hand, in the gap emerge now QP solutions 
with very weak Zp <C 1 which reconnect the FS into a 
large one. We are dealing nevertheless with effectively 
truncated FS with well developed arcs. The effect of 
larger k is essentially to increase Zp in the gaped region, 
in particular near (jr, 0). Finally, for large K = 0.6 which 
corresponds to the regime consistent with optimal doping 
or overdoping in cuprates, Zp is essentially only weakly 
decreasing towards (tt, 0) and the FS is well pronounced 
and concave as naturally expected for t' < 0. 




mm 



max 



FIG. 1. (color) Contour plot of spectral functions 
j4(k,oj = 0) at T — for various k in one quarter of the 
Brillouin zone. 

In order to understand the pseudogap features at low 
but finite k we present in Figs. 2,3 A(k,ui) for k = 0.1. 
Spectra in Fig. 3 are presented along the lines a - d in 
the Brillouin zone as shown in Fig. 2. As expected from 
Eq. (H^) the pseudogap is smallest along the zone diag- 
onal (line a) where, moreover, the pseudogap appears at 



oj > 0, so that it would not be seen in ARPES. Lines 
a and b are thus examples of the region where arcs of 
the FS are well pronounced, i.e., their QP weight is not 
strongly renormalized, Zp < Z ci . On the other hand, 
following lines c and d the chemical potential u> = falls 
into the pseudogap. We see in Fig. 3(a), (b) that QP in 
fact crosses coherently the FS {w = 0) although with very 
small Zp C 1. It is evident from Fig. 1 that Zp within 
the pseudogap remains small only for k <C 1 while it 
increases and finally smears out the concept of the pseu- 
dogap for k > k* ~ 0.5. 



(0#) 




(0,0) (mM 

FIG. 2. (color) A(k,w = 0) for n = 0.1 and lines a) - d) 
used in Fig. 3. 




-0.4 0.0 0.4 0.8 -0.8 -0.4 0.0 0.4 

FIG. 3. A(k,uj) for k — 0.1 along different directions in the 
Brillouin zone, corresponding to Fig. 2. 
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It is quite remarkable to notice that in spite of Zp <C 1 
the QP velocity vf is not diminished within the pseudo- 
gap. In fact it is even enhanced, as seen in Fig. 3(c) 
and even more clearly in Fig. 4 where the contour plot of 
A(k, uj) is shown corresponding (with restricted k span) 
to Fig. 3(c). Again, it is well evident in Fig. 4 that QP is 
well defined at the FS, while it becomes fuzzy at uj ^ 
merging with the solutions E^, respectively, away from 
the FS. 




(7r/8,57r/8) 



(3;r/8,77r/8) 



FIG. 4. (color) Contour plot of A(k,uj) for k = 0.1, cor- 
responding to Fig. 3(c). Note large quasiparticle velocity 
around uj ~ 0. 

Presented formalism offers a possible scenario for the 
evolution of the FS with doping from a pocket-like into 
a large one. In order to explain results in Figs. (3), (4) 
on the effective truncation of the FS and the character 
of QP within the pseudogap we note that it is essentially 
enough that both k and uj k are finite to yield well defined 
FS. Since gross features do not depend on the particular 
form Eq. ( p9| ) we present here simplified analysis using 
instead 



X "(Q + q,^) = 



C[S(uj - uj k ) - 5(u + w K )],q_L < k, 
0, q± > k, 



(38) 



where q±_ denotes the component perpendicular to the 
AFM zone boundary. Let us assume that e£ f - fi = e ~ 
/i = e ~ 0. We also linearize dispersion e£. f 
(1,1), so that we get from 



and 



k-Q 



at the FS and take that v£ f 



Eq. (£51), 



A 2 (w + UJ K + 6 - Uj)(uJ K + Uj) 

2w (w + uj k — e + uj) (uj k — uj) 



(39) 



where w = v^K and A = A^. Let us evaluate QP prop- 
erties on the FS assuming that it is located at e = 0, i.e., 
on the AFM zone boundary. We obtain from Eq. ( |39| ) 
the QP weight Zp 



Z cf 



■ 9S' 




-1 








duj 


w=0,e=0- 





uj k (uj k + w) 



(40) 



This clearly leads to and explains Zp <C 1 for uj k w ~ 
2v^Jn 2 <§; A 2 . This is generally the case within the 
gaped part of the FS for small k < k* , as shown in 
Fig. 1. It should be also noted that the latter condition 
is essentially always satisfied near (ir, 0) where v£ f ~ 
and consequently also w ~ 0. 

Let us evaluate in the same way the QP renormalized 
velocity vf at the FS. Here we realize that the k depen- 
dence of £' is essential. The latter is in Eq. ([l(j]) given 
by the e dependence, 



VF 

as' 



(i 



dsy z F 
!k>~zrt : 



A 2 



uj=0,e=0 CJ K ((jJ K + w) ' 



(41) 



which in contrast to Zp leads to an enhancement of vf- 
In the case uj k w <C A 2 we thus get 



VF 



U>h 
W 



2J 



(42) 



Final vf is therefore not strongly renormalized, since 2 J 
and are of similar order. Furthermore, vp is enhanced 
in the parts of FS where is small, in particular near 
(tt, 0) point. The situation is thus very different from 
'local' theories where E(k,w) ~ £(w) and the QP renor- 
malization is governed only by Zp. In our case the 'non- 
local' character of S(k, uj) is essential in order to properly 
describe QP within the pseudogap region. 

Let us discuss further the behavior of the DOS jV(uj), 
Eq. (|3C|). It is evident from Fig. 1 that the contribu- 
tion to jV(uj ~ 0) will come mostly from FS arcs near the 
zone diagonal while the gaped regions near (tt, 0) will con- 
tribute less. Results in Fig. 5 (full lines) the development 
of jV(uj) with k, as corresponding to FS in Fig. 1. Wc 
see that the DOS indeed reveal a pseudogap at uj < A 
however the pseudogap is visible only for k < 0.5 and 
deepens for k — > 0. 



S 0.4 



().() 
0.6 
S 0.4 
0.2 
0.0 



- K=0.04 A 

NX 1 s 


- K=0.2 ^ 






1 1 1 1 1 1 1 

K=0.4 


■ 1 1 1 1 1 ■ 

K=0.6 

^ \ 
-C \ 

r / ^-»«^_ s 
/ / \ ~~ 




/ 

r / — 



-0.4 



0.4 



2 0.0 0.2 -0.4 -0.2 0.0 C 

co/r o)// 
FIG. 5. Density of states jV(uj) (full lines) and weighted 
DOS jV w (lu) (dashed lines) for different k. 
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The DOS is measured in cuprates via angle integrated 
PES, e.g. for LSCO in flpj , as well as via the scanning 
tunneling microscopy (STM) . It is well possible that 
within both experiments the matrix elements are essen- 
tial leading to enhanced contribution near the x points 
in the Brillouin zone. It has been proposed that for the 
c-axis conductivity ]36| the interplanar hopping should 
be weighted by the matrix element 



tu(k) 



(cos fc T — cos k„ 



(43) 



The same arguments as for the c-axis conductivity might 
apply also for the STM effective DOS as well as for inte- 
grated PES, therefore we present also the weighted DOS 
N w where w(k) is introduced additionally into Eq. (30). 
Results also presented in Fig. 5 (dashed line) show much 
stronger pronounced pseudogap, in particular at low k. 
This is quite evident since u>(k) essentially destroys the 
effect of FS arcs near (tt/2,tt/2) which present the main 
contribution (due to small velocity in hole-pocket FS) to 
the usual N(uj). 

In Fig. 6(a) we show the average Z av along the FS, as 
well as the QP DOS, defined as 



aZ 



cf 



QP 



2tt 2 



cISf 



(44) 



in Fig. 6(b) we present as well as the dependence of DOS 
at the FS, both N(0) and N w (0), as a function of n. 
Note that Nqp should be relevant for the specific heat, 
i.e., Nqp oc 7 = CV/T at low T (provided that we are 
dealing with a normal Fermi liquid). It is quite impor- 
tant to understand that decreasing k (smaller doping) 
also means the decrease of Nq p , which is consistent with 
the observation of the pseudogap also in the specific heat 
in cuprates J37j. We note here that such a behavior is 
not evident when one discusses the metal-insulator tran- 
sition. Namely, in a Fermi liquid with (nearly) constant 
Fermi surface one can drive the metal-insulator transition 
by Z^v — > and within an assumption of a local charac- 
ter Ti(oj) this would lead to vf —>■ and consequently to 
Nqp — + 00. Clearly, the essential difference in our case 
is that within the pseudogap regime £(k,u;) is nonlocal, 
allowing for the simultaneous decrease of N(0) and Nqp- 
Finally, let us comment on the influence of finite T. 
While T enters within this approach also via effective 
parameters as u£ f and predominantly k(T), here we con- 
sider only the direct effect via the thermodynamic factor 
in Eq. (p5|). It is evident that T > smears out Eif. 
This becomes important at small k in particular for QP 
in the pseudogap regime. In Fig. 7 we present j4(k, cj), 
corresponding to Fig. 3(c), for several values of low T. 
The main conclusion is, that weak (but sharp) QP peak 
with Zp < 1 at T = is smeared out already by very 
small T > T s ~ 0.02 t and is not at all visible (there is 
no overdamped peak) at higher T, the remainder being 
an incoherent background at ui ~ for T > T s . This is 



important to realize that ARPES experiments in fact do 
not observe no well defined QP peak near (ir, 0) in the 
underdoped regime at T > T c . 
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FIG. 6. (a) Average QP weight Z av and QP DOS Nqp vs. 
k. (b) DOS N(0) and weighted DOS AC(0) vs. k. 
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-0.4 0.0 0.4 -0.4 0.0 0.4 -0.4 0.0 0.4 

00/f CQ/f C0/f 

FIG. 7. A(k,u) for k = 0.1 and k along the cen- 
tral part of the line c in Fig. 2 for various T: T = 0, 
T = 0.02t, and T = 0.05*. Momentum k ranges from 
(7r/4,37r/4) =p (7r/32,7r/32). 
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V. SELF CONSISTENT CALCULATION 

The full selfconsistent set of equations for E = E pm + 
En, Eqs.(H),(||), and for G, Eq.©, is solved numer- 
ically. For given /j, the FS emerges as a solution deter- 
mined by the relation Ck F + E'(ki?,0) = fi. We should 
note that at given /i, electron concentration c e as cal- 
culated from the DOS Af(uj), Eq. (|3l|) , does not in gen- 
eral coincide with the one evaluated from the FS vol- 
ume, c e = Vfs/Vq. Nevertheless, apart from the fact 
that within the t-J model validity of the Luttinger the- 
orem is anyhow under question psj], in the regimes of 
large FS both quantities appear to be quite close. The 
position of the FS is mainly determined by Ck and E pm , 
while in this respect Ey is less crucial. 

As discussed in Sec. Ill we use in Eq. (Q) the most ap- 
propriate and simplest approximation to insert the un- 
renormalized ^4°(k, u>), i.e., the spectral function without 
a self-consistent consideration of Eif but with E pm fully 
taken into account. We choose here t' = —0.2t and again 
k = \fcfi while 771 and r\i are determined as a function 
of Ch from model calculations We use TV = 40 x 40 
points in the Brillouin zone and broadening e/t = 0.05. 

In Fig. 8 we present hole concentration Ch vs. /1 as 
obtained from Af(cu) at T = 0. We solve selfconsistent 
equations by iteration, whereby for 0.06 < < 0.11 we 
find in the equations an instability signaled by oscilla- 
tory behavior instead of the convergence and an unique 
solution can not be obtained in the region indicated by 
the dashed line. However, at lower (and higher) doping 
the solution is converged. It seems that the region of in- 
stability coincides with the transition from the large to a 
small FS. 



i(k) = 



A(k, uj)duj. 



(45) 




FIG. 8. Hole doping Ch as a function of the chemical po- 
tential n/t, following from the selfconsistent calculation. 

The shape of the FS is most clearly presented with con- 
tour plot of the electron momentum distribution function 
defined as 



Results for a characteristic development of the FS with 
Ch are shown in Fig. 9. At higher doping, Ch = 0.26 
and also Ch — 0.20, we get a common large FS topol- 
ogy. In the intermediate doping regime, ch = 0.14, the 
pseudogap is pronounced at momenta around (71", 0) and 
the FS shows a tendency of forming a small FS. The gap 
is more pronounced because of longer AFM correlation 
length £ (smaller k). At Ch < ~ 0.06 solutions are 
consistent with a small pocket-like FS whereby this be- 
havior is enhanced by t' < as realized in other model 
studies H . On increasing doping the FS rather abruptly 
changes from a small into a large one as suggested from 
the results of the SCBA |39|. The smallness of c° has 
the origin in quite weak dispersion dominated by J and 
if at Ch — * which is overshadowed by much larger £k at 
moderate doping, where the FS is large and its shape is 
controlled by t'/t. 



(Ojr) 





(0.0) (/r,0)(0,0) (/r,0) 

FIG. 9. (color) Electron momentum distribution n(k) for 
various ch- Thin contour lines represent n(k) in increments 
0.1 while heavy line is corresponding to h(k) = 0.8 for 
c h = 0.26, 0.20, 0.14 and n(k) = 0.85 for c h = 0.04. 

In Fig. 10(a) and Fig. 10(b) we present calculated 
A(k, uj) along the principal directions in the Brillouin 
zone, i.e., (tt/2,tt/2) ~ * (0,0) — > (tt,0). It is evident that 
E pm leads to a strong damping of hole QP and quite in- 
coherent momentum-independent spectrum A(k, uj) for 
w <C — J which qualitatively reproduces ARPES and nu- 
merical results Jl2|| . Electron QP (at u > 0) are in gen- 
eral very different, i.e., with much weaker damping aris- 
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ing only from S pm . Note a relatively high QP velocity in 
the higher doping regime Ch = 0.26, Fig. 10(a), as com- 
pared to a more narrow dispersion on the scale 2 J at low 
doping Ch — 0.04 , Fig. 10(b), where we find the regime 
of small pocket-like FS. 



S 



(a) J=03t 
f'=-0.2f 
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-3 -2 
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FIG. 10. A(k, uj) along main directions in the Brillouin 
zone. 

In Fig. 11 we present the development of spectral func- 
tions at fixed Ch = 0.2, but now varying k as an indepen- 
dent parameter. Let us concentrate on the emergence of 
the pseudogap near (ir, 0). At n = 0.4 ~ yfch the pseu- 
dogap is essentially not yet developed. Nevertheless, the 
gap opens with decreasing k, in particular for (at this 
doping unrealistic value) k — 0.05. 
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0.3 
(fift 

-0,3 
-0.6 



0.6 
0.3 
Gift 

-0.3 
-0.6 







K=0.4 





- k=02 



Mm 



(0,0) 




(/r,0) (0.0) 
] max 

FIG. 11. (color) Contour plot of A(k, u>) for fixed doping 
Ch = 0.2 and various k. 

VI. CONCLUSIONS 

We have presented the theory of spectral functions 
within the t-J model whereby our method is based on 
the EQM for projected fermionic operators and on the 
decoupling approximation for the self energy assuming 
the fermions and spin fluctuations as essential coupled 
degrees of freedom. We first make some comments on 
the method: 

a) The EQM approach for spectral functions (as well 
as for other dynamical quantities) seem to be promising 
since it can treat exactly the constraint which is essential 
for the physics of strongly correlated electrons. 

b) In finding the proper approximation for the self energy 
within the EQM approach it is plausible that the main 
coupling is of fermions with spin fluctuations, where close 
to the AFM ordered state both transverse and longitu- 
dinal spin fluctuations are important. It is well possible, 
however, that other contributions could be important, 
e.g., the coupling to pairing fluctuations. 

c) In an ordered AFM our method reproduces naturally 
the results for the spectral function (of a hole) within 
SCBA, which is highly nontrivial, since both approaches 
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are quite different. 

d) The coupling to longitudinal spin fluctuations appear 
to be most important for QP near the AFM zone bound- 
ary and is responsible for the opening of the pseudogap. 
Here the coupling is only moderately strong and can be 
treated in the lowest-order decoupling scheme. 

e) The present theory uses the spin response as an input. 
x"(q, (J), Eq. (p9|) , corresponds in general to a Fermi liq- 
uid or to a short-range AFM liquid. Results remain in 
fact qualitatively similar as far as x' '(<!> ^-0 is nonsingular. 
In the opposite case, e.g., if we would use the MFL form 
as an input, the Fermi surface would still be defined, but 
QP would have vanishing weight Z — > [p2| . 

Let us further discuss some main results of the pre- 
sented theory: 

a) The fermion-paramagnon coupling as manifested in 
S pm remains effective and strong even at moderate dop- 
ing. The full calculation shows that the coupling leads 
to a large incoherent part in the hole part [u> < 0) of the 
spectral function, as well as to the renormalized band e£ f • 

b) The main consequence of due to the coupling to 
longitudinal spin fluctuations is the appearance of the 
pseudogap at k < k*. The pseudogap opens predomi- 
nantly along the AFM zone boundary and its extent is 
qualitatively given by Eq.(|37|), dependent on J and t' but 
not directly on t. Evidently the pseudogap has similarity 
to d- wave-like dependence along the FS, for t' < being 
largest near the (it, 0) point. 

c) How strong is the pseudogap effect depends mainly 
on k. At small k <C k* parts of the Fermi surface near 
(ir/2, 7r/2) remain well pronounced (for t' < 0) while the 
Fermi surface within the pseudogap is suppressed, i.e., 
QP have a small weight Zp <C 1, in particular near zone 
corners (w, 0). 

d) The simplified analysis yields large Fermi surface, al- 
though a truncated one, except at very small k <C k* 
where Sit by itself induces a small hole-pocket-like Fermi 
surface. On the other hand, E pm generates hole pock- 
ets already for Ch < c° ~ 0.06. In fact, an instability 
of the self consistent calculation indicates the emergence 
of hole-pocket Fermi surface even at c/, < 0.1. However, 
it is well possible that within the present approximation 
scheme E pm is overestimated at intermediate c/j, an in- 
dication for that being quite a weak dispersion . 

e) Our method is approximative in the evaluation of S, 
hence it is not surprising that the volume of the Fermi 
surface does in general not coincide with the one follow- 
ing the Luttinger theorem. It is anyhow questionable if 
such a relation should be valid within the t-J model |3|] 
due to the projected basis and strong correlations. Nev- 
ertheless, in the regime of a large Fermi surface the full 
calculation yields the Fermi surface volume quite close to 
the Luttinger one. 

f) For k < k* the QP within the pseudogap have small 
weight Zp <C 1 but not diminished vp, which is the effect 
of the nonlocal character of S(k, w). A consequence is 



that QP within the pseudogap contribute much less to 
QP DOS Nqp- This can explain the reduction of the 
latter with doping and the appearance of the pseudogap 
in the specific heat being essential for the understanding 
of the specific heat in underdoped cuprates. 
g) Although most results are presented for T = 0, we 
can discuss some effects of T > 0. First effect is that 
within the pseudogap the QP with Zp -C 1 are washed 
out (not just overdamped) already for very small T < 
T s <C A PG . On the other hand, the pseudogap is mainly 
affected by k. So we can argue that the pseudogap would 
be observable for «(c/i,T) < k* ~ 0.5. This effectively 
determines the pseudogap crossover temperature T*(ch). 
From the quantitative studies of the t-J model it 
appears that in the region of interest k is nearly linear in 
both T and Ch so we would get approximately 

T* ~ To (1 — Ch/ c£), (46) 

where T * - 0.6J and c* h ~ 0.15. 

Finally we make some comments on the relevance of 
our results to experiments on cuprates, in particular with 
respect to observed pseudogap and Fermi surface fea- 
tures: 

a) The (large) pseudogap scale shows in ARPES on 
BSCCO as a hump at ~ lOOeV §. Our results indi- 
cate quite a similar pseudogap scale, e.g. in Fig. 5 the 
u) < pseudogap ~ 0.3 t (note that t ~ 0.4eV), since 
A PG is determined mainly by J and t'. 

b) The truncated Fermi surface in underdoped BSCCO 
appears as an arc (part of the large Fermi surface corre- 
sponding to t' < 0) in the Brillouin zone ||, effectively 
not crossing the AFM zone boundary, which is also char- 
acteristic of our results for k < K*, originating from the 
strong coupling to spin fluctuations with commensurate 
(jr, 7r). The same is the origin of shadow features in spec- 
tral functions pronounced at intermediate doping and in 
particular at weak doping. 

c) Our results for the depletion of the DOS N(0) (N w (0)) 
with decreasing doping are qualitatively consistent with 
the integrated PES (so far known for LSCO |M) and 
STM J3j|, although in this relation the importance of 
matrix element corrections is not yet clarified. The 
same holds for the calculated decrease of QP DOS Nqp 
with doping essential in connection with the specific- 
heat pseudogap in underdoped cuprates |37|) . It should 
be however mentioned that our results for both Af(0) 
as well as of Nqp oc 7 indicate on weaker suppression 
with decreasing doping than observed in experiments. 
This is due to remaining contribution of Fermi surface 
arcs, which could be overestimated in our approach for 

k < K*. 

d) Both the value and the dependence of the pseudo- 
gap temperature T*(c^), as estimated in Eq.(^), seem 
to be very reasonable in connection with experimental 
evidence, arising from various transport and magnetic 
properties in cuprates W. 
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